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A 610-MHz survey of the Lockman Hole with the Giant Metrewave 
Radio Telescope - 1. Observations, data reduction and source 
catalogue for the central 5 deg 2 
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1 INTRODUCTION 



The Lockman Hole dLockman. Jahoda & McCammort 1986) is the 
region of lowest Hi column density in t he sky, with the low i nfrared 
background (0.38 MJy sr -1 at 100 urm lLonsdale et alj2003h mak- 
ing this region particularly well suite d for deep infrare d observa- 
tions. The Spitzer Space Telescope dWerner et ai]|2004) observed 
~12 deg 2 of the region in 2004 as par t of the Spitzer Wide -area In- 
frared Extragalactic survey (SWIRE: lLonsdale et alj|2003l) . Obser- 
vations were centered on 10 h 45 m 00 s , +58°00'00" (J2000 coordi- 
nates, which are us ed throughout this work) using the Infrared Ar- 
ray Camera (IRAC; lFazio et all2004T) operating at 3.6, 4.5, 5.8 and 
8 jlm, and the M ultiband Imaging Photometer for Spitzer (MIPS; 
iReike et alj2004h at 24, 70 and 160 fim. 

A great deal of complementary data has been taken on the 
Lockman Hole at other wavelengths in order to exploit the avail- 
ability of sensitive infrared observations. There have been deep 
optical observations in the U, g', r' and i 1 bands taken with the 
Mosaic-I camera at Kitt Peak National Observatory (KPNO), to 
5<T Ve ga magnitude limit s of 24.1, 25.1, 24.4 and 23.7 respec- 
tively dSurace et alj 120051) . There is existing near-infrared data 
across the region from the Two Micron All Sky Survey (2MASS; 
iBeichman etai1l2003l) . to J, H and A" s band magnitudes of 17.8, 
16.5 and 16.0, and a band-merged catalogue containing 323 044 
sources from the IR AC, MIPS 24 jj.m 2MASS and KPNO data 
has been produced dSurace et alj I2005T) . A photometric redshift 
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Lockman Hole has been con structed from the band-merged data 
dRowan-Robinson et al . 2008). Data Release Six of the Sloan Digi- 
tal Sky Survey (SPSS DR6; lAdelman-McCartv et alj2007h covers 
the whole region in the ugriz bands, with both photometric and 
spectroscopic observations available to a depth of ~22 mag. Fur- 
ther deep infrared observations of the Lockman Hole are under- 
way as part of the Deep Extragalact ic Survey section of t he UK 
Infrared Deep Sky Survey (UKIDSS; [Lawrence et aT]|2007l) in the 
J, H and K bands, with a target sensitivity of J ~ 22 mag. There 
have been deep surveys of the Lockman Hole wit h the Submil- 
limet re Common- User Bolometer Arr ay (SCUBA; iHolland et al.1 
Il999h at 850 urn dCoppin et alj|2006h . and with t he X-ray satel- 
lites ROSAT dHasinger et al j|l998h. XMM-Newton dHasinger et al.1 
200 ll ; Mainieri et al. 2002); iBrunner et alj 1200 8h and Chandra 



(Polletta etal. 20C 



A variety of radio surveys cover areas within the Lockman 
Hole. The first of these was by Ide Ruiter et al. I dl997h , who ob- 
served 149 1.4-GHz sources within an area of 0.35 deg 2 , using 
the Very Large Array (VLA) in C-configuration, with an rms noise 
level o f 30-55 [iJy beam - . A similar deep observation was carried 
out by [Ciliegi et al.1 d2003l) . who observed 63 sources at 4.9 GHz 
within a 0.087 deg 2 region, using the VLA in C-configuration, 
with an rms noise l evel of 11 /ijy beam -1 . More recently, 
Biggs & Ivisoi] (2006) found 506 sources within a 0.35 deg area, 
using the VLA at 1 .4 GHz operating in the A- and B-configurations, 



catalogue containing 229 238 galaxies and quasars within the and with an rms noise level of 4.6 /ijy beam -1 . The Faint Images 

of the Radio Sky at Twenty-cm (FIRST; iBecker. White & Helfand 
1 19951) and NRAO VLA Sky Survey (NVSS; ICondon et aljfl998l) 

* E-mail: tsg25@cam.ac.uk surveys both cover the entire region at 1.4 GHz, but only to rel- 
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Right Ascension (J30O0) 



Figure 1. The arrangement of 12 pointings observed in the central region 
of the Lockman Hole. 




Figure 2. The uv coverage for pointing 10. Baselines less than 1 kA were 
not used in imaging and have been omitted from the figure. 



tively. 

There is a clear need for deep radio observations over a 
significant area within the Lockman Hole, in order to extend 
the multi-wavelength information presently available for a large 
number of sources. In this paper, we present observations of 
the Lockman Hole taken at 61 MHz with the Gian t Metre- 
wave Radio Telescope (GMRT; lAnanthakrishnaiJ 1 20051) . cover- 
ing ~5 deg 2 of sky with a resolution of 6 x 5 arcsec 2 , posi- 
tion angle (PA) +45° and centred on 10 h 45 m 00 s , +58°00'00" 
to match the SWIRE coverage. This is the third in a series 
of 610-MHz GMRT surveys targeting the Spitzer deep legacy 
survey regions, following the S pitzer extragalactic First Look 
Survey (xFLS: iGarn et al] 12007) and the Euro pean Large Are a 
ISO Survey-North 1 (ELAIS-N1) survey (ENh lGarn etalll2008h . 
This work, in combination with the deep optical and infrared 
data, will be used to study the infrared / radio correlation 
for st a r-form ing galaxie s (e.g. Helou, Soifer & Rowan-Robinson 
19851: ICondonf Il992l; iGarrettl 120021: lAppleton et al.J |2004|; 



Murphy et all 120061: IVlahakis Eales & Dunnd 120071 : 
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20081 : IGarn. Ford & Alexandej 120081) . and the link between star 
formation and Active Galactic Nuc lei (AGN) activity (e.g. 
iRichards et al]|2007l : iBundv et alj|2007l) . as wel l as the pro perties 
of the faint radio source population at 610 M Hz iBond fet al .120071 : 
iMoss et all2007l : iTasse et alj2007l : IGarn et alj|2008h . 

In Section [2] we describe the observations and data reduction 
techniques used in the creation of the survey. Section[3]presents the 
mosaic and a source catalogue containing 2845 sources above 6(7. 
In Section|4]we calculate the 610-MHz differential source counts, 
and using local radio luminosity functions obtained from the liter- 
ature, we demonstrate that a model of pure luminosity evolution is 
sufficient to account for the differential source counts at 610 MHz 
and 1 .4 GHz simultaneously. 

Throughout this work a flat cosmology with £1^ = 0.7 and 
Hq = 70 km s Mpc~' is assumed. 



2 OBSERVATIONS AND DATA REDUCTION 

The 'central' region of the Lockman Hole, consisting of 12 point- 
ings spaced by 36 arcmin in a hexagonal pattern (shown in Fig. [TJ 
was observed in 2004 July 24 and 25, using the GMRT operating 
at 610 MHz. A further 26 pointings surrounding the central region 
(the 'outer' region) were observed in 2006 July in order to complete 
the coverage of the majority of the region observed by Spitzer. The 
outer observations suffered from a series of power failures in one 
of the arms of the GMRT, and required several catch-up observing 
sessions to cover the whole field. This paper describes the obser- 
vations of the central 12 pointings only, with the results from the 
outer region to be presented in Paper II. 

The flux density scale was set through observations of 3C48 
or 3C286 at the beginning and end of each observing session. The 
task SETJY was used to calculate 610-MHz flux densities of 29.4 
and 21.1 Jy respectively, using the Astro nomical Image Pro cessing 
Software (AIPS) implementation of the iBaars et al.l ( 119771) scale. 
Each field was observed for a series of interleaved scans, varying 
between 8 and 10 min in duration, in order to maximise the uv cov- 
erage. The scans took place over a total time period of 7 and 9 hours 
on the two days, with the typical total time being spent on each field 
being 48 min. A nearby phase calibrator, J1035+564, was observed 
for 4 min between every three scans to monitor any time-dependent 
phase and amplitude fluctuations of the telescope. The measured 
phase typically varied by less than 10° between phase calibrator 
observations. 

Observations were made in two 16-MHz sidebands centred 
on 610 MHz, each split into 128 spectral channels, with a 16.9 s 
integration time. An error in the timestamps o f the uv data wa s 
corrected using UVFXT (for further details see IGarn et ai] |2007t) . 
Initial editing of the data was performed separately on each side- 
band with standard AIPS tasks, to remove bad baselines, antennas, 
and channels that were suffering from large amounts of narrow 
band interference, along with the first and last integration periods 
of each scan. The flux calibrators were used to create a bandpass 
correction for each antenna. In order to create a continuum chan- 
nel, five central frequency channels were combined together, and an 
antenna-based phase and amplitude calibration created using obser- 
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vations of J1035+564. This calibration was applied to the original 
data, which was then compressed into 1 1 channels of bandwidth 
1.25 MHz, each small enough that bandwidth smearing was not a 
problem for our images. The first and last few spectral channels, 
which tended to be the noisiest, were omitted from the data, lead- 
ing to an total effective bandwidth of 13.75 MHz in each sideband. 
Further flagging was performed on the 1 1 -channe l data set, and the 
two sidebands combined using UVFLP (again, see lGarn et al .l2007h 
to improve the uv coverage. The uv coverage for one of the point- 
ings is shown in Fig. [2] Baselines shorter than 1 kX were omitted 
from the imaging, since the GMRT has a large number of small 
baselines which would otherwise dominate the beam shape. 

Each pointing was divided into 31 smaller facets, arranged in 
a hexagonal grid, which were imaged individually with a separate 
assumed phase centre. The large total area imaged (a diameter of 
l.°8, compared with the full width at half-maximum of the GMRT, 
which is ~0.°74) allows bright sources well away from the observed 
phase centre to be cleaned from the images, while the faceting pro- 
cedure avoids phase errors being introduced due to the non-planar 
nature of the sky. All images were made with an elliptical restoring 
beam of size 6x5 arcsec 2 , PA +45°, with a pixel size of 1.5 arcsec 
to ensure that the beam was well oversampled. 

The images went through three iterations of phase self- 
calibration at 10, 3 and 1 min intervals, and then a final iteration of 
phase and amplitude self-calibration, at 10 min intervals, with the 
overall amplitude gain held constant in order not to alter the flux 
density of sources. The self-calibration steps improved the noise 
level by about 10 per cent, and significantly reduced the residual 
sidelobes around the brighter sources. 

The final rms noise before correction for the GMRT primary 
beam was between 57 and 63 ilJy beam -1 for pointings 1 to 9 
and pointing 12, with higher noise levels of around 80 tlJy beam -1 
for the remaining two pointings. The expected thermal noise limit 
isj~55 /iJy beam -1 , calculated using Equation 1 from lGarn et al.l 
d2007l) . The increased noise level for pointings 10 and 1 1 is due to 
a very bright radio source (FIRST J103333. 2+581502), with 1.4- 
GHz flux density of ~2 Jy, located ~0?75 from the pointing cen- 
tres. This source lies within the region cleaned during the imaging 
process, but residual sidelobes remain in our mosaic, and the in- 
creased noise level is visible in Fig. [3] 

In our prev ious 610-MHz GMRT surveys jGarn et alj|2007l : 
iGarn et we detected a systematic difference in the mea- 

sured flux densities of sources present in more than one point- 
ing. This is thought to be due to an elevation-dependent error in 
the position of the GMRT primary beam (N. Kantharia, privat e 
communication). We repeated the analysis of IGarn et al.l d2007h . 
correcting each point ing for the nominal primary beam shape 
dKanfharia & Raoll200ll) and modelling the effective phase centre 
of the primary beam as having a systematic offset from its nominal 
value. A shift of ~2.5 arcmin in a north-east direction was required 
for each pointing in order to remove the systematic flux density 
offset. 

The 12 pointings were mosaicked together, taking into account 
the offset primary beam and weighting the final mosaic appropri- 
ately by the relative noise of each pointing. The mosaic was cut 
off at the point where the primary beam correction dropped to 
20 per cent of its central value, a radius of 0P53 from the centre 
of the outer pointings. The final mosaicked image is available via 
http://www.mrao.cam.ac.uk/surveys/, and a sample region of the 
mosaic, away from bright sources, is shown in Fig. [4] to illustrate 
the quality of the image. 




10 54 52 50 48 46 44 42 40 38 36 
RIGHT ASCENSION (J2000) 



Figure 3. The rms noise of the final mosaic, calculated using S EXTRACTOR 
(see SectionfJJ. The grey-scale ranges between 60 and 300 (I Jy beam - 1 , and 
the contours are at 100 and 200 ti Jy beam -1 respectively. The increased 
noise for pointings 10 and 1 1 (top right and centre-right) is due to a very 
bright radio source located at 10 h 33 m 33.2 s , +58°15'02". 



3 SOURCE CATALOGUE 

Source Extractor (SEXTRACTOR; iBertin & Arnoutsl 1 19961) was 
used to calculate the rms noise a across the mosaic. A grid of 
16 x 16 pixels was used in order to track changes in the local 
noise level, which varies significantly near the brightest sources. 
Fig. [3] illustrates the local noise, with the grey-scale varying be- 
tween 60 /iJy beam -1 (the approximate noise level in the centre of 
the pointings) and 300 /iJy beam -1 (the expected noise level at the 
distance where the GMRT primary beam gain was 20 per cent of its 
central value). The 100 and 200 /iJy beam -1 contours are plotted 
on Fig. [3] to guide the eye. Our GMRT data suffer from dynamic 
range issues near the brightest sources, and the final mosaic has 
increased noise and residual sidelobes in these regions, as for the 
xFLS and EN1 surveys. While the local noise calculated by SEX- 
TRACTOR increases due to these residuals, some of them still have 
an apparent signal-to-noise level that is greater than 6(7. We there- 
fore opted for a two-stage selection process for our final catalogue. 

3.1 Catalogue creation 

An initial catalogue of 5223 sources was created using SEXTRAC- 
TOR. The mosaic was cut off at the point where the primary beam 
correction dropped to 20 per cent of its central value, but only 
sources inside the 30 per cent region were included in the cata- 
logue to avoid the mosaic edges from affecting the local noise esti- 
mation. The requirements for a source to be included in the initial 
catalogue were that it had at least 5 connected pixels with bright- 
ness greater than 3 times the local noise a, and a peak brightness 
greater than 6o". The image pixel size meant that the beam was rea- 
sonably oversampled, so the source peak was taken to be the value 
of the brightest pixel within a source. The integrated flux density 
was calculated using the FLUX.AUTO option within SEXTRAC- 
TOR, which creates an elliptical aperture around each object (as 
described in lKronlll98d) . and integrates the flux contained within 
the ellipse. 

Artefacts are seen near bright sources which may be in- 
cluded in the initial catalogue. We used the technique described 
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Figure 4. A 0.2 deg 2 region of the 610-MHz mosaic, illustrating the image quality. The grey-scale ranges between —0.2 and 1.0 mjy beam and the noise 
is relatively uniform and around 70 jUly beam -1 , apart from small regions near bright sources where the noise increases. The brightest source in the field is 
GMRTLH J 105043.2+575231, with a peak brightness of 50.4 mjy beam -1 . 



in iGarn et all d2008h in order to identify regions that are affected 
by these artefacts, and within these areas required sources to have a 
peak brightness greater than 12(7 in order to be included in the final 
catalogue. All sources with peaks greater than 10 mjy beam -1 were 
examined for artefacts. The final catalogue contains 2845 sources - 
we have erred on the side of caution in order to produce a catalogue 
with little contamination from spurious sources. 

Table [T] presents a sample of 10 entries in the catalogue, 
which is sorted by RA. The full table is available as Supple- 
mentary Material to the online version of this article, and via 
http://www.mrao.cam.ac.uk/surveys/ Column 1 gives the IAU des- 
ignation of the source, in the form GMRTLH Jhhmmss.s+ddmmss, 
where J represents J2000.0 coordinates, hhmmss.s represents RA in 
hours, minutes and truncated tenths of seconds, and ddmmss rep- 
resents the Dec. in degrees, arcminutes and truncated arcseconds. 
Columns 2 and 3 give the RA and Dec. of the source, calculated 
using first moments of the relevant pixel brightnesses to give a cen- 
troid position. Column 4 gives the brightness of the peak pixel in 
each source, in mjy beam -1 , and column 5 gives the local rms 
noise in /xjy beam -1 . Columns 6 and 7 give the integrated flux 
density and error, calculated from the local noise level and source 
size. Columns 8 and 9 give the X, Y pixel coordinates from the mo- 
saic image of the source centroid. Column 10 is the S EXTRACTOR 
deblended object flag - 1: where a nearby bright source may be af- 
fecting the calculated flux, 2: where a source has been deblended 
into two or more components from a single initial island of flux, 
and 3: when both of the above criteria apply. There are 77 sources 
present in our catalogue with non-zero deblend flags, and it is nec- 
essary to examine the images to distinguish between the case where 
one extended object has been represented by two or more entries, 
and where two astronomically distinct objects are present. 
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Figure 5. Radio spectral index a between 610 MHz and 1.4 GHz, for 
sources in the FIRST catalogue. 



3.2 Comparison with the FIRST survey 

In order to test the positional accur acy of our catalogu e, we paired 
it with sources in the FIRST survey teecker et al.ll995h . The whole 
of our 610-MHz survey region is covered by FIRST, and 340 
sources are found with positions within 6 arcsec in the two surveys. 
The difference in source positions in the GMRT catalogue relative 
to the VLA FIRST survey is approximately Gaussian, with mean 
offset of 0.01 arcsec in RA and —0.1 arcsec in Dec. The standard 
deviations of the distribution are 0.3 and 0.5 arcsec respectively. 
We make no alteration to the coordinates of sources in our survey 
due to the good agreement in positions. 

The spectral index distribution of the matched sources is 
shown in Fig. [5] using the integrated flux density measurements. 
The distribution peaks around a = 0.8, where a is defined so that 
the flux density 5 scales with frequency v as 5 = 5()V -0! . There are 
significant biases in this distribution due to the differing sensitivity 
levels of the two surveys - Fig. [6] shows the variation in spectral 
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Table 1. A sample of 10 entries from the 610-MHz Lockman Hole catalogue, sorted by ascending RA. The full version of this table is available as Supple- 
mentary Material to the online version of this article, and via http://www.mrao.cam.ac.uk/surveys/ 



Name RA Dec. Peak Local Noise Int. Flux Density Error X Y Flags 

J2000.0 J2000.0 mjybeam" 1 ^iJybeam"' mJy mJy 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


GMRTLH J104237.0+583759 


10:42:37.04 


+58:37:59.1 


0.564 


82 


1.666 


0.136 


3994 


4773 





GMRTLH J104237. 1+572634 


10:42:37.10 


+57:26:34.7 


0.539 


88 


0.473 


0.072 


4018 


1917 





GMRTLH J104237.4+574652 


10:42:37.41 


+57:46:52.1 


0.445 


68 


0.412 


0.056 


4010 


2729 


3 


GMRTLH J104238.3+572630 


10:42:38.36 


+57:26:30.7 


0.698 


88 


0.586 


0.082 


4012 


1914 





GMRTLH J104239.2+585011 


10:42:39.27 


+58:50:11.5 


0.536 


73 


0.626 


0.085 


3978 


5261 





GMRTLH J104239.3+572752 


10:42:39.35 


+57:27:52.1 


0.545 


78 


0.643 


0.081 


4006 


1969 





GMRTLH J104240.6+580353 


10:42:40.62 


+58:03:53.6 


0.438 


71 


0.786 


0.083 


3987 


3409 





GMRTLH J104240.6+573836 


10:42:40.65 


+57:38:36.5 


0.410 


60 


0.331 


0.054 


3995 


2398 





GMRTLH J 10424 1.4+5 80505 


10:42:41.46 


+58:05:05.9 


0.648 


77 


0.710 


0.082 


3982 


3458 





GMRTLH J104241.6+573259 


10:42:41.63 


+57:32:59.6 


0.547 


78 


0.518 


0.078 


3992 


2173 






s 
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Figure 6. The variation in spectral index a with 610-MHz flux density. 
Crosses represent sources in the GMRT and FIRST catalogues, with the 
solid black line showing the limiting spectral index that could be detected, 
given the respective sensitivity levels. The 610-MHz flux density limit is 
shown by the dotted black line. 



index with the 610-MHz flux density. Below a few mJy, steep spec- 
trum sources that are visible at 610 MHz cannot be detected in the 
FIRST survey, due to its limited depth (shown by the solid line). 



3.3 Sub-millimetre galaxies 

Bigas & Ivison (2007) have produced a sample of 12 sub- 
millimetre galaxies (SMGs) detected with combined VLA and 
Multi-Element Radio-Linked Interferometer (MERLIN) observa- 
tions at 1.4 GHz. The majority of sources have a flux density at 

1.4 GHz less than 150 /iJy, but there are three 'bright' sources, 
SMG06, SMG10 and SMG11 with 1.4-GHz flux densities of 246, 
295 and 249 /xJy respectively, and typical errors of ~ 10 /iJy. One 
of these (SMG10, at z = 1.212) is detected in the 610-MHz cata- 
logue with a flux density of 538 ± 83 jl Jy, giving it a spectral index 
of a = 0.72 ±0.19. The remaining two radio-bright SMGs, and all 
of the faint ones, are undetected at 610 MHz. The 610-MHz noise 
level at the location of the two undetected bright SMGs is approx- 



imately 100 /xjy beam , leading to an upper limit on the spectral 
index of these sources of a = 1.05. 



4 DIFFERENTIAL SOURCE COUNTS 

We calculated source counts for the Lockman Ho le using a metho d 
similar to the one used for the ELAIS-N1 field iGarn et alj|2008h , 
by deriving faint and bright source counts from two separate re- 
gions of the mosaic. In order to calculate source counts below 
10 mJy, we chose a 75 x 75 arcmin 2 region that was relatively 
free from bright sources above 10 mJy (and consequently from 
the errors resulting from those sources). For reliability, source 
counts were constructed from all sources with a signal-to-noise ra- 
tio greater than 7, and binned by their integrated fl ux density, using 
the same flux density bins as in iGarn et alj l lioolh . The lowest flux 
bin went down to 556 flly, corresponding to the approximate noise 
of 80 fl Jy beam - 1 within the selected region. We calculated the dif- 
ferential source count for sources with flux densities above 10 mJy 
by considering the central ~3.3 deg 2 of the map, excluding the re- 
gion near pointing 10 which had significantly greater noise. Within 
this area, residual sidelobes are seen from several bright sources, 
but inspection of the residuals found them to have a typical flux 
density of below 1 mJy, which will not affect the source counts. 

The differential source counts dN /dS were corrected for the 
fraction of the image over which they could be detected, taking 
into account the increase in noise near the bright sources. In or- 
der to correct for resolution bias, we applie d a correction fac tor of 
3 per cent for source counts below 1 mJy (Moss et al. 20Q3). No 
correction for resolution bias was made for brighter sources. Ta- 
ble [2] gives the source counts, mean flux density (5) of sources in 
each bin, d/V/dS and d/V/dS normalised by (S) 2 ' 5 , the value ex- 
pected from a static Euclidean universe, along with Poisson error 
estimates calculated from the number of sources detected in each 
flux density bin. The source counts are consistent with our results 
from t he Spitzer extragalactic First Look Survey and ELAIS-N1 
fields iGarn et alj|20071 ; IGarn et alj|2008h . and in Table [3] we give 
the source counts found by combining data from the three fields. 

Fig. |7] shows the normalised differential source counts from 
our three survey fields, alo ng with 610-MHz source counts from 
two other GMRT surveys (Mo ss et alj|2007l ; iBondi et alj|2007h . 
and a collection of shallow survey s made with the Westerbork 
Synthesis Radio Telescope (WSRT; IValentiin et afll 19771 ; iKatgertl 
1 19791 : IValentiinlll980l : [Katgert -Merkeliin et al.lll985f) . The turnover 
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Table 2. 610-MHz differential source counts for the Lockman Hole survey. 



r.2.5 



Flux Bin 
(mJy) 



(S) N 
(mJy) 



dN/dS 



dN/dS (S) 2 



(sr~' Jy -1 ) (sH Jy LS ) 



0.556- 


0.798 


0.674 


94 


108.8 


9.4±1.0x 10 8 


1 1 .1 ± 1.2 


0.798 - 


1.145 


0.951 


61 


61.3 


3.7±0.5 x 10 8 


10.4±1.3 


1.145- 


1.643 


1.381 


31 


31.0 


1.3±0.2x 10 8 


9.3± 1.7 


1.643 - 


2.358 


1.992 


27 


27.0 


7.9 ± 1.5 x 10 7 


14.0 ±2.7 


2.358- 


3.384 


2.909 


26 


26.0 


5.3 ± 1.0 x 10 7 


24.3 ±4.8 


3.384- 


4.856 


3.981 


14 


14.0 


2.0±0.5 x 10 7 


20.0 ±5.3 


4.856 - 


6.968 


5.700 


9 


9.0 


8.9 ± 3.0 x 10 6 


21.9±7.3 


6.968 - 


10.00 


8.429 


6 


6.0 


4.2 ± 1.7 x 10 6 


27.1 ±11.1 


10.00- 


13.49 


11.70 


14 


14.0 


4.0± 1.1 x 10 6 


59.3± 15.9 


13.49- 


18.20 


15.70 


8 


8.0 


1.7±0.6x 10 6 


52.5±18.6 


18.20- 


24.56 


20.71 


17 


17.0 


2.7 ± 0.7 x 10 6 


165 ±40.0 


24.56 - 


33.14 


27.65 


11 


11.0 


1.3 ±0.4 x 10 6 


163 ±49.2 


33.14- 


44.72 


37.70 


7 


7.0 


6.0 ±2.3 x 10 s 


167±63.1 


44.72 - 


60.34 


57.47 


2 


2.0 


1.3±0.9x 10 5 


101 ±71.7 


60.34 - 


81.41 


74.23 


4 


4.0 


1.9 ± 1.0 x 10 5 


285 ±142 


81.41 - 


109.8 


94.89 


4 


4.0 


1.4±0.7x 10 s 


390 ±195 


109.8 - 


148.2 


120.5 


3 


3.0 


7.8±4.5 x 10 4 


395 ±228 


148.2- 


200.0 


160.3 


1 


1.0 


1.9±1.9x 10 4 


199 ±199 



Table 3. Average 610-MHz differential source counts for the Spitzer extra- 
galactic First Look, ELAIS-N1 and Lockman Hole surveys. 



AN 
dS 



.-1.5 



n <h(s,z)h(z)(i 

Jo 



zfdz 



(1) 



Flux Bin 


(S) 


N 


N c 


dN/dS 


dN/dS (S) 2 ' 5 


(mJy) 


(mJy) 






(sr- 1 Jy -1 ) 


(sr" 1 Jy L5 ) 


0.270 - 


0.387 


0.331 


269 


399.6 


3.8±0.2x 10 9 


7.6±0.5 


0.387 - 


0.556 


0.463 


269 


287.2 


1.9±0.1 x 10 9 


8.8±0.5 


0.556 - 


0.798 


0.669 


268 


288.3 


8.7±0.5 x 10 8 


10.1 ±0.6 


0.798 - 


1.145 


0.958 


172 


172.3 


3.6±0.3 x 10 8 


10.3 ±0.8 


1.145- 


1.643 


1.365 


125 


125.0 


1.8±0.2x 10 8 


12.6 ±1.1 


1.643- 


2.358 


1.951 


84 


84.0 


8.6±0.9x 10 7 


14.5±1.6 


2.358- 


3.384 


2.793 


73 


73.0 


5.2±0.6x 10 7 


21.5±2.5 


3.384- 


4.856 


4.060 


53 


53.0 


2.6 ± 0.4 x 10 7 


27.7 ±3.8 


4.856 - 


6.968 


5.712 


33 


33.0 


1.1 ±0.2 x 10 7 


28.2±4.9 


6.968 - 


10.00 


8.414 


28 


28.0 


6.8±1.3x 10 6 


43.9±8.3 


10.00- 


13.49 


11.55 


56 


56.0 


5.5 ±0.7 x 10 6 


78.6±10.5 


13.49- 


18.20 


15.77 


38 


38.0 


2.8±0.5 x 10 6 


86.1 ±14.0 


18.20- 


24.56 


21.01 


49 


49.0 


2.6 ± 0.4 x 10 6 


168 ±24.1 


24.56 - 


33.14 


28.57 


36 


36.0 


1.4±0.2x 10 6 


198±33.0 


33.14- 


44.72 


38.16 


24 


24.0 


7.1 ± 1.4 x 10 s 


202 ±41. 2 


44.72 - 


60.34 


54.58 


9 


9.0 


2.0 ± 0.7 x 10 5 


137±45.7 


60.34- 


81.41 


68.26 


11 


11.0 


1.8±0.5 x 10 5 


217±65.5 


81.41 - 


109.8 


94.89 


4 


4.0 


4.8 ± 2.4 x 10 4 


133 ±66.7 


109.8 - 


148.2 


121.6 


6 


6.0 


5.3 ± 2.2 x 10 4 


276 ±113 


148.2- 


200.0 


175.1 


4 


4.0 


2.6±1.3x 10 4 


339 ±170 



where <j) z (L) has been separated into luminosity evolution <j>o(S,z) 
and density evolution (1 +z) F terms, and the change in comoving 
volume w ith redshift is given by h(z) - for more details, see Sec- 
tion 4.3 of lSevmour et all d2004b . 

The three populations of radio sources consist of a 'steep' 
spectrum AGN population with assumed spectral index of a = 0.8, 
a 'flat' spectrum AGN population with spectral index of and a 
'starburst' population, also assu med to have spectral index of 0.8. 
Following ISevmour et all d2004l) we consider no density evolution 
(P = 0), and use the 1.4 -GHz AGN luminosity functio ns and lu- 
minosity evolution from iRowan-Robi nson et all d 19931), which i s 
based on earlier work at 2.7 GHz by iDunlop & Peacock! (l990). 
The AGN contribution to the radio source counts is therefore fixed 
at the same values found by previous authors, and dominates the 
source counts above a few mJy. 

We obtain a local RLF for starburst galaxies from 

iMauch & Sadleil d2007l) . measured from 4006 local radio sources 
present in the 6dF Galaxy Survey ( I Jones et all |2005|) and NVSS 
surveys, class ified as starbur s t gala xies through their optical spec- 
tra. Following lHopkins et al.1 d 19981) the starburst luminosity evolu- 
tion is parameterised as (1 ±z)*s where we assume evolution takes 
place between z = and 2, based on measurements of the change in 
cosm i c star formation rate across this redshift range (Mad au et al.1 
1 1996) ; iHopkins & Beacoml 120060 . The star formation rate is as- 
sumed to remain constant for z > 2. An M82-type galaxy at z = 2 
would have a flux density below 1 fi Jy at 1 .4 GHz, and would there- 
fore not contribute to the observed source counts, making this as- 
sumption valid. 

We construct models for the source counts at 610 MHz and 
1.4 GHz, with the 1.4-GHz RLFs converted to 610 MHz through 
the use of the spectral indices chosen above. Fig. [7] shows the pre- 
dicted 610-MHz source counts from three models of pure lumi- 
nosity evolution, with the starburst evolution parameter Q taking 
the values 2, 2.5 and 3, in order to compare to previous works. The 
source coun t measurements at 1.4 GHz are taken f r om a wide range 
of surveys dWhite et al.||l997l; Iciliegi et alj|l999l; Gruppioni et al. 



1999; Prandoni et al. 2001; Bondi et al. 2003; Hopkins et al. 2003; 
Seymour et all |2004| ; iHuvnh et all 120051 ; iBiggs & Ivisod 12006); 
Simps on et alj|20060 . and are shown in Fig. [8] along with the same 



three luminosity evolution models. 

We find that a value of Q between 2.5 and 3 agrees well 
with th e data. From 1 ,4 -GHz source counts, [R owan-Robinso n et al.1 
dl993l) found O = 3.1.ICondon. Cotton & Broderickl d2002l) found 
Q = 3 ± 1, [Seymour et alj d2004h found O = 2.5 ± 0.5 and 
Huvn h et al.l d20050 foundO = 2.7. Using 1.4-GHz inform ation 



in source counts is clearly visible below ~2 mJy, d ue to the emer- 
gence of a new population of radio sources (see e.g. lSimpson et al.1 
120060 . 



and a range of other star formatio n indicators, [Hopkins (2004]) find 
Q = 2.7 ± 0.6 and P = 0. 15 ± Q.6. IMqss et al.1 d2007l) use their 610- 
MHz source counts to fit a value of Q = 2.451Q 25 . Our results are 
therefore consistent with previous findings, and demonstrate that 
both sets of source counts can be modelled simultaneously with the 
same radio luminosity functions and evolutionary dependence. 



4.1 Source count models 

Radio source counts have been use d to understand the evolution of 
radio sources for many years (e.g. iLongaial 19660 . We model the 
source counts as b eing due to three distinct populatio ns of radio 
sources, following Seymour, McHardv & Gunn (2004). The nor- 
malised differential source counts for a population with radio lu- 
minosity function (RLF) given by (j) z (L) can be calculated using 



5 CONCLUSIONS 

We have observed ~5 deg 2 of the Lockman Hole with the GMRT 
at 610 MHz. The majority of the observations have a rms noise of 
~60 /xJy beam -1 before primary beam correction, with the noise 
level increasing in the direction of a bright FIRST source outside of 
our image. A catalogue of 2845 radio source components detected 
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1 1 ■ — ■ — ■ — ■ — 1 — ■ — ■ — 1 — 1 

10" 4 O.OOl 0.01 0.1 1 

Flux (Jy) 

Figure 7. Differential source counts at 610 MHz iValentiin et alj|l977t iKatgerll [l979l ; IValentiinl Il98d; iKatgert-Merkeliin et alj|l985t iBondi et al.l |2007| ; 
iMoss et al ] |2007l ; lGarn etai]|2007l ; iGarn et alj [2008). normalised by the value expected in a static Euclidean universe. Model source counts calculated with 
pure luminosity evolution with Q = 2 (dotted line), 2.5 (dashed line) and 3 (solid line) are shown - see Section |4~Tl for more details. 



- x White (1997) 




10" 5 10" 4 0.001 0.01 0.1 1 

Flux (Jy) 



Figure 8. Different i al source counts at 1.4 GHz Jwhite et alj|l997l; ICiliegi et alj|l999l; [Gruppioni et alj 1 19991; IPrandoni et at] 1200 ll ; IBondi et alj |2003| ; 
Hopkins et al. 2003; Seymour et al. 2004; Huvnh et al. 2005; Biggs & Ivison 2006; Simpson et al. 2006), normalised by the value expected in a static Eu- 
clidean universe. Model source counts calculated with pure luminosity evolution with Q = 2 (dotted line), 2.5 (dashed line) and 3 (solid line) are shown - see 
Section l4~T1 for more details. 



above 6<J has been presented, and the full catalogue and mosaic 
made available via http://www.mrao.cam.ac.uk/surveys/ 

We calculate differential source counts at 610 MHz by con- 
sidering two separate regions of the mosaic. The source counts are 
calculated between 556 fl]y and 200 mJy, and are shown to agree 
with previous 610-MHz source counts in other regions. We present 
average source counts calculated from this work and two previ- 
ous surveys of the Spitzer extragalactic First Look Survey field 
and ELAIS-N1 field. The turnover in differential source counts is 
clearly visible below 2 mJy, and we show that a three-component 
population containing steep spectrum AGN, flat spectrum AGN and 



starburst galaxies which undergo pure luminosity evolution is suf- 
ficient to model the 610-MHz and 1.4-GHz source counts simulta- 
neously. 
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